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Nucleon Decay Matrix Elements for Domain- Wall Fermions 
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We report on the nucleon decay matrix elements with domain-wall fermions in quenched approximation. Results 
from direct and indirect method are compared with a focus on the process of a proton decaying to a pion and a 
lepton. We discuss the renormalization necessary for the matching to the continuum theory. Preliminary results 
for the renormalized chiral lagrangian parameters are presented. 



1. INTRODUCTION 

Nucleon decay is one of the most important 
aspect that any (SUSY) GUT model has. In 
low-energy effective theories it is represented 
by dimension six operators made of one lepton 
field and three quark fields. In the decay of 
p — > 7T° + e + , for example, we need to calculate 
hadronic matrix element of the three quark oper- 
ator C| /L;L = e i]k {u' T CP R/L di)P L u k with ini- 
tial proton and final pion states, where C is the 
charge conjugation matrix. Pr/l represents the 
right or left projection matrix, respectively. 

Lattice calculation of the nucleon decay matrix 
elements historically started with the indirect cal- 
culation (I). First, one estimates the overlap of 
the operator to the proton state, a and (3, 

aP L u p = (0|of.>;fc = 0>, (1) 
f3P L u p = (0|of ;L |p;fc = 0), (2) 

with Up being the proton spinor. Then, with the 
help of chiral perturbation theory, one can calcu- 
late the nucleon decay matrix elements. The first 
attempt to calculate the nucleon decay matrix ele- 
ment directly handling the three point function is 
done by Gavela et. al ||, showing significant de- 
viation from the indirect calculations. Recently 
JLQCD || pointed out the incompleteness of the 
calculation The matrix element has a tensor 
structure, 

^P\°l/L; L \P^) = Pl[W - iiW q ]u p , (3) 
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where q = k — p. The relevant form factor Wq 
is what we need since the term proportional to ^ 
vanishes after multiplying the lepton spinor. W q 
is called the irrelevant form factor. JLQCD uses 
quenched Wilson gauge configurations and Wil- 
son fermion at parameters summarized in Table ^ 
(c). The three quark operators are renormalized 
with one-loop perturbation theory. Now the dif- 
ference of direct and indirect calculations is not 
huge, but, still 30 to 40% for most cases. 

In this study we try to calculate nucleon decay 
matrix elements with domain- wall fermions in the 
quenched approximation. Using DBW2 gauge ac- 
tion makes chiral symmetry breaking especially 
small compared to the other actions With 
this good chiral symmetry one expects the good 
property to calculate hadronic matrix elements: 
1) preventing operator mixing with different chi- 
ral structure, and 2) good scaling, even down to 
1/a ~ 1 GeV region ||||. We restrict ourselves 
to the case of degenerate quark mass in the me- 
son. In this case, one can calculate p — > tt° + l + 
and p — > 7r + + v decay amplitudes. The latter is 
obtained by multiplying the former by \[2 under 
the exact SU(2) symmetry of u and d quarks. 



2. PROCEDURE AND RESULTS 

The parameters of our simulation are shown in 
Table |l| (a). Our quark propagator is obtained 
by averaging two propagators calculated using 
periodic and anti-periodic boundary conditions. 
Thus, the effective temporal extent is 64. 100 in- 
dependent gauge configurations are analyzed for 
the three and two point functions to calculate the 
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Table 1 

Lattice parameters and corresponding inverse lattice spacings with p mass input. Values of ro are listed 
for another estimate of the scale. We use the parameter set (a) for the nucleon decay, f and \ represent 
unpublished data with 50 and 100 configurations respectively. 



fcrmion size (V x T) L s AJ 5 l/a p [GeV] reT gauge 6/y r /a 



ref. 



16 3 " x 32 
16 3 x 32 



(a) DWF 

(b) DWF 



12 

1G 



1.8 

1.8 



28 3 x 48 x i 
16 3 x 32 



1.23(5) 
1.31(4) 



± 



DBW2 

DBW2 



0.87 

0.87 



3.58(4) { 



(c) Wilson 

(d) DWF 



16 1. 



2.30(4) 
1.92(4) 



Wilson 
Wilson 



6.0 
6.0 



5.368 



ratio 



(Mti)Q% /L;L (t)J P (t )) 



(J*(ti)Jl(t)){J p (t)J p (to) 



(4) 



In the three point function proton and pion in- 
terpolating fields are located at to = 6 and 
t\ = 24 respectively. Momentum ±p with pa = 
(1,0, 0)7r/8 or (1, 1,0)tt/8 is injected to the pion 
and the operator in the three point function, as 
well as in the pion two point function in the de- 
1 Z-x and \fZ~n are overlap of 



nominator. 

and J p to the corresponding pion and proton 
states, which is estimated from the fit of two point 
functions. We take four different quark masses 
m f a = 0.02, 0.04, 0.06, 0.08, among which the 
strange mass is located between the last two. Af- 
ter appropriate projection, trace and the subtrac- 
tion of the irrelevant contribution, Wq is calcu- 
lated with a fit to plateau of the ratio. The re- 
suits with bare operator C^. L are shown in Fig- 
ure [I]. Wo can only depend on q 2 and on mj 
which can enter through the masses of proton and 
pion. We assume a form of the fitting function of 
Wo = Co + ciq 2 + c 2 (q 2 ) 2 + c 3 m/. In our preci- 
sion the physical kinematics are approximated as 
rrif — and — q 2 = 0. The filled diamond in the 
figure represents the result extrapolated to this 
point. We also calculate the chiral lagrangian pa- 
rameters a and j3 extrapolated to m/ =0. By 
the tree-level chiral perturbation theory ||, the 
relevant part of the matrix element is obtained, 

(7T°|0f.» rel ~ 0(1 +D + F)P L u p /V2f, (5) 

where D — 0.47 and F — 0.80 from the exper- 
iment H, and / is f w = 0.131 GeV. The open 
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Figure 1. Wo for p — > (ir°,l + ) decay with Ol-l 
as a function of invariant mass squared of lepton. 

diamond in Fig. [j] represents this indirect estima- 
tion. There is no significant difference between di- 
rect and indirect results in the current statistics. 
However, 40% excess for the central value from 
indirect method is quite similar to the JLQCD 
result EH. Similar of result is obtained also for 

(Ao% L \ P ). 

The above results need be renormalized for the 
value in the MS scheme in continuum theory. We 
expect these operators are multiplicatively renor- 
malized as the chiral symmetry breaking is small 
(m res ~ 1.5 MeV). The renormalization factor 
for DBW2 gauge action is available for one-loop 
perturbation theory fic| with L s treated as in- 
finity. In Table @ wc list the values of renormal- 
ization factors for the nucleon decay by pertur- 
bation theory with two different mean field (MF) 
improvement schemes. The MF factor u is either 
P 1 / 4 or 1/8K C , where P stands for plaquette, K c 



3 



Table 2 

Rcnormalization factors for the axial vector cur- 
rent and the nucleon decay three-quark operator 
in MS with NDR. Scale p = 1/a is chosen for 
nucleon decay. 



operator 


perturbative JlfJ 
plaq K c 


non 
perturbative 


W L;R/L 


0.976 0.769 
1.009 0.713 


0.7776 (5) 



the critical hopping parameter of the four dimen- 
sional Wilson fermion. We also list the renor- 
malization factor for the axial vector current, for 
which we know the non-perturbative value at the 
parameter point of (b) in Table [l] with the method 
using the conserved current || . For the axial vec- 
tor current renormalization, the mean field im- 
provement with K c is in better agreement with 
the non-perturbative one, while the other is 30% 
larger. This difference is considered as the sys- 
tematic error of the perturbative estimate. The 
difference arises mainly from the quark wave func- 
tion renormalization and the MF factor. Thus, 
this problem does not apply to quantities defined 
by a particular ratio in which those factors cancel 
out. Bk is one of those cases. 

In Table | we list the result of a and (3 of ours 
and by JLQCD. The operators are renormalizcd 
at p, = 1/a with MS for both results. Renormal- 
ization factor with K c is employed for our calcu- 
lation. Our values of a and (3 in [GeV 3 ] obtained 
by setting the scale with the p mass input, are 
quite different from those with Wilson fermion 
by JLQCD. The origin of the discrepancy could 
be, 1) discretization error, 2) systematic error of 
the MF perturbation, 3) difference of the renor- 
malization scale. 2) is about 30% for DWF, 3) 

Table 3 

Chiral lagrangian parameter renormalized at p = 
1/a in MS in NDR. Scale from p mass is used for 
results in [GeV 3 ]. Quoted errors are statistical 
only. 





[GeV^] 


x 10" 2 








\a\ 




|ot|rg 


wi 


This work 


0.6(1) 


0.7(1) 


0.15(3) 


0.18(4) 


JLQCD@ 


1.5(1) 


1.4(1) 


0.19(2) 


0.18(2) 



is almost negligible. 1) is expected to be most 
severe as it is seen in the difference of the scale 
from p mass for Wilson and DWF with Wilson 
gauge background at same 1/g 2 (see (c) and (d) 
in Table |l|). Indeed, if we set the scale from dif- 
ferent quantity , which eliminates order a error 
in the scale for Wilson fermion, the values are 
consistent. 

3. SUMMARY 

We calculated proton decay matrix elements 
with DBW2 gauge action and DWF. We did not 
find a significant difference in values obtained 
by direct and indirect methods within our cur- 
rent numerical precision. Further simulation with 
non-degenerate quark masses will reduce the er- 
ror of direct calculation and make the extrac- 
tion of nucleon to kaon decay amplitudes possi- 
ble. Within the current systematic error of the 
perturbation theory the chiral lagrangian param- 
eters for the nucleon decay calculated with DWF 
are consistent with those obtained with Wilson 
fermion by JLQCD if the scale is set by ro . More 
stringent comparison will be done for those val- 
ues after non-perturbative renormalization is per- 
formed. 
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